function plot_outputs(LineNumber,model_name)
switch model_name
    case 'human_mean'
        load('/home/shu/workspace/Research/2DKnee_Torso_Outputs/data/IndexAll.mat')
        % choose load which matrix
        str_linenumber = num2str(LineNumber);
        amat_name = [str_linenumber 'a_opt.mat'];
        
        load(['/home/shu/workspace/Research/2DKnee_Torso_Outputs/result/' amat_name]);
        a0 = a_opt;
        % choose the output combination
        hiptype = IndexAll(LineNumber,1);
        slopetype = IndexAll(LineNumber,2);
        torsotype = IndexAll(LineNumber,3);
        load('HM_1014.mat');
        ns_time = mean(time_rec)*ave_timeSK';
        
    case 'NAO'
        load('/home/shu/workspace/Research/NAO_Output_Test/data/IndexAll.mat')
        % choose load which matrix
        str_linenumber = num2str(LineNumber);
        amat_name = [str_linenumber 'a_opt.mat'];
        
        load(['/home/shu/workspace/Research/NAO_Output_Test/result/' amat_name]);
        a0 = a_opt;
        % choose the output combination
        hiptype = IndexAll(LineNumber,1);
        slopetype = IndexAll(LineNumber,2);
        torsotype = IndexAll(LineNumber,3);
        load('NAO_1014.mat');
        ns_time = mean(time_rec)*ave_timeSK';
        
    case 'AMBER'
        load('/home/shu/workspace/Research/AMBER_Output_Test/data/IndexAll.mat')
        % choose load which matrix
        str_linenumber = num2str(LineNumber);
        amat_name = [str_linenumber 'a_opt.mat'];
  
        load(['/home/shu/workspace/Research/AMBER_Output_Test/result/' amat_name]);
        a0 = a_opt;
        % choose the output combination
        hiptype = IndexAll(LineNumber,1);
        slopetype = IndexAll(LineNumber,2);
        torsotype = IndexAll(LineNumber,3);
        load('AMBER_1014.mat');
        ns_time = mean(time_rec)*ave_timeSK';
    case 'fitNAO'
        load('/home/shu/Documents/Research_Related/NAO_CODE/Data/Afit.mat')
        % choose the output combination
        hiptype = 2;
        slopetype = 2;
        torsotype = 1;
        load('NAO_1014.mat');
        ns_time = mean(time_rec)*ave_timeSK';
end

% for hip related output
if hiptype == 1
    meanValuePlotG1(ave_timeHP,upperBHP,lowerBHP,x_meanHP,'HipPos.eps',...
        'Position(m)',a0(1,:),1,ns_time,1,LineNumber);
elseif hiptype == 2
    meanValuePlotG1(ave_timeHPL,upperBHPL,lowerBHPL,x_meanHPL,'HipPosLinear.eps',...
        'Position(m)',a0(1,:),1,ns_time,2,LineNumber);
end


% for non-stance slope related outputs
if slopetype == 1 % non-stance slope
    meanValuePlotG1(ave_timeNSslope,upperBNSL,lowerBNSL,x_meanNSL,'nsSlope.eps',...
        'Slope',a0(2,:),2,ns_time,3,LineNumber)
elseif slopetype == 2 % linearized non-stance slope
    meanValuePlotG1(ave_timeNSslopeL,upperBLNS,lowerBLNS,x_meanLNS,'nsSlopeLiear.eps',...
        'Slope',a0(2,:),2,ns_time,4,LineNumber)
elseif slopetype == 3 % hip angel
    meanValuePlotG1(ave_timeHip,upperBHip,lowerBHip,x_meanHip,'HipAngle.eps',...
        'Angle(rad)',a0(2,:),2,ns_time,5,LineNumber);
end

% for torsorelated outputs
if torsotype == 1 % torso hip angle
    meanValuePlotG1(ave_timeTH,upperBTH,lowerBTH,x_meanTH,'TorsoAngle.eps',...
        'Angle(rad)',a0(5,:),2,ns_time,8,LineNumber);
elseif torsotype == 2 % stance COM slope
    meanValuePlotG1(ave_timeSCOM,upperBSCOM,lowerBSCOM,x_meanSCOM,'stCOM.eps',...
        'Slope',a0(5,:),2,ns_time,9,LineNumber);
elseif torsotype == 3 % linearized stance COM slope 
        meanValuePlotG1(ave_timeSCOML,upperBSCOML,lowerBSCOML,x_meanSCOML,'stCOMLinear.eps',...
        'Slope',a0(5,:),2,ns_time,10,LineNumber);
% % % elseif torsotype == 2 % theta 4
% % %     meanValuePlotG1(ave_timeT4,upperBT4,lowerBT4,x_meanT4,'Theta4.eps',...
% % %         'Angle(rad)',a0(5,:),2,ns_time,9,LineNumber);
% % % elseif torsotype == 3 % non-stance torso slope
% % %     meanValuePlotG1(ave_timeNST,upperBNST,lowerBNST,x_meanNST,'NSTorsoSlope.eps',...
% % %         'Slope',a0(5,:),2,ns_time,10,LineNumber);
% % % elseif torsotype == 4 % linearized non-stance torso slope
% % %     meanValuePlotG1(ave_timeNSTL,upperBLNST,lowerBLNST,x_meanLNST,'NSTorsoSlopeLinear.eps',...
% % %         'Slope',a0(5,:),2,ns_time,11,LineNumber);
% % % elseif torsotype == 5 % stance torso slope
% % %     meanValuePlotG1(ave_timeST,upperBST,lowerBST,x_meanST,'STorsoSlope.eps',...
% % %         'Slope',a0(5,:),2,ns_time,12,LineNumber);
% % % elseif torsotype == 6 % linearized stance torso slope
% % %     meanValuePlotG1(ave_timeST,upperBLST,lowerBLST,x_meanLST,'STorsoSlopeLinear.eps',...
% % %         'Slope',a0(5,:),2,ns_time,13,LineNumber);
% % % elseif torsotype == 7 % com
% % %     meanValuePlotG1(ave_timeCOM,upperBCOM,lowerBCOM,x_meanCOM,'COM.eps',...
% % %         'Position(m)',a0(5,:),2,ns_time,14,LineNumber);
% % % elseif torsotype == 8 % linearized com
% % %     meanValuePlotG1(ave_timeCOML,upperBCOML,lowerBCOML,x_meanCOML,'COMLinear.eps',...
% % %         'Position(m)',a0(5,:),2,ns_time,15,LineNumber);
end

% stance knee
meanValuePlotG1(ave_timeSK,upperBSK,lowerBSK,x_meanSK,'SKnee.eps',...
    'Angle(rad)',a0(3,:),2,ns_time,6,LineNumber);
% non stance knee
meanValuePlotG1(ave_timeNSK,upperBNSK,lowerBNSK,x_meanNSK,'NSKnee.eps',...
    'Angle(rad)',a0(4,:),2,ns_time,7,LineNumber)

% mkdir([model_name 'figure']);
movefile('figure/*.eps',[model_name 'figure'])
close all;
end